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Electron correlation in cuprates leads to a global constraint 
Ak = on the gap function Ak resulting in a gap nodal 
surface. We give physical arguments supported by numerical 
results and discuss some experimental results to argue that 
correlations also lead to a local constraint on charge fluctua- 
tions in k-space close to the Fermi surface, which may result 
in a substantial overlap of the Fermi surface with the gap 
nodal surface. 



An actively debated issue in the experimental study of 
high-T c cuprate superconductors p], p3| , p!6[ is the symme- 
try and detailed k-dependence of the gap function Ak. 
We suggest a form of Ak, based on physical arguments 
and numerical results, which is in qualitative agreement 
with the available experimental results. The form that 
we suggest is a generalised version of an early suggestion 
of Anderson || in the context of cuprates and even an 
earlier suggestion of Cohen || in the context of sodium 
metal. In our scenerio, what is of primary importance 
is the vanishing of gap function on patches of the fermi 
surface(FS), leading to a gapless superconducting state 
with finite density of states for quasiparticle excitations 
at the chemical potential. That is, in the vicinity of finite 
portions of the FS the gap function has a non-zero value 
Ak ~ Ai and in the vicinity of the rest of the portions on 
the FS the gap function has the form Ak ~ (|k| — ki?)A 2 . 
Here Ai and A 2 are finite and are weakly dependent on 
the wave vector k. The issue of the symmetry of the 
gap function in the context of the above suggestions is 
discussed at the end of this paper. 

The local constraint on double occupancy in the con- 
ducting CuC>2 layer of the cuprates has non-trivial conse- 
quence on the structure of the gap function. In a simple 
tight binding one band large-[/ Hubbard model or the 
t — J model of the CuOi layers (for which there is am- 
ple theoretical and experimental evidence) this implies 
essentially zero on-site pairing amplitude in the super- 
conducting state: A^ = {c\^c\^) — 0, which in k-space 
means 



assuming usual zero momentum pair condensation. The 
above is a global constraint, which forces the gap function 



to develop a nodal line (in 2D) or nodal surface (in 3D), 
such that the integral of Ak on the two sides of the nodal 
surface cancel each other. We call this a gap nodal surface 
(GNS). 

An important question is the shape and location of the 
gap nodal surface in k-space. A simple possibility is a gap 
with d x 2_ y 2 symmetry g which has nodal lines k x = ±k y 
and the global constraint (Eq. (1)) is ensured by the 
symmetry. On the other hand, the originial RVB scenerio 
of Anderson || as well as a simple mean field theory || 
for the undoped Mott insulator suggested a gap function 
(for the pre-existing neutral Cooper pairs) of the form 
cos fcr+cos k y = 0, which satisfies the global constraint in 
k-space and also the GNS coincides with the tight binding 
FS at half filling. Anderson has suggested that on doping 
the Mott insulator the GNS will also shrink and continue 
to overlap with the FS, thereby making it a truly gapless 
superconductor. However, Anderson did not give any 
physical arguments in support of his conjecture 

In this letter, we will argue that real space correlation 
changes the nature of charge fluctuations in k-space close 
to the FS and it is incompatible with the nature of charge 
fluctuations demanded by a finite value of Ak on the FS. 
This incompatibility reduces the pairing correlations (if 
possible, to zero) on the FS. We also present numerical 
results on small one and two dimensional Hubbard and 
t — J clusters to substantiate the form of our gap function. 
Towards the end, we discuss some experimental results 
including the recent ARPES results in the light of 
our suggestion. 

Let us begin by examining the behaviour of the stan- 
dard BCS wave function close to the FS: 
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where the product is over only half of the k-space (e.g. 
k x > and all k y in 2D), with & t _ k = -^{c^c_ k ^ — 

cj^clj^) as the singlet pair creation operator on points 
(1) k and — k in k-space. It is clear that u k is the proba- 
bility amplitude of finding zero singlet pair with momen- 
tum (k, — k), v^ttk^k is that of finding one singlet pair 
(of charge 2e), and v\ that of finding two singlet pairs 
(of total charge 4e). The BCS state has identical phase 
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relations for various configurations of pair occupancy in 
k-space. That is, when the product in Eq. (|3[) is ex- 
panded out, the resulting sum has identical phase for all 
terms, each term corresponding to different configura- 
tions of the (k, — k) occupancy. Superconductivity can 
thus be thought of as a coherent charge-2e fluctuating 
state in k-space. Since itk^k is non-zero only in a thin 
energy shell around the FS, the coherent 2e charge fluc- 
tuation is concentrated around the FS. (It is interesting 
to note that this coherence in k-space results in phase 
coherence among the Cooper pairs in real space also). 
Away from the shell, we either have a completely filled 
band (inside the FS) or a completely empty band (out- 
side the FS) and hence no charge fluctuations. Thus, 
MkWk(~ Ak) is a measure of coherent charge fluctuations 
or k-space electron pair compressibility. 

We now argue that strong correlations in real space 
lead to a suppression of such coherent charge fluctuations 
close to the FS. Strongly correlated electrons in ID and 
2D, described by a large-C/ Hubbard model, have certain 
unique features close to the FS. It is well known that in 
the ID Hubbard model there is singular forward scat- 
tering between two electrons with opposite spins close 
to the FS. This leads to a finite phase shift || at the 
FS and the consequent failure of the Fermi liquid the- 
ory, resulting in the vanishing of the discontinuity in rik 
at the FS (Luttinger liquid behaviour). It also implies 
an effective hard-core repulsive pseudopotential between 
electrons with opposite spins close to the FS. Thus no two 
electrons close to the FS, with opposite spins, can have 
the same momentum, thereby making k-points close to 
the FS essentially singly occupied. Single occupancy in 
the vicinity of the FS has been seen in numerical works 
on one and two dimensional t — J models |?]]. Ours is a 
first and natural explanation of this. Freezing of occu- 
pancy to one at every point in k-space close to the FS 
also means reduced (or vanishing) pair fluctuations on 
the FS. Coherent pair fluctuations are thus unlikely to 
develop on or very close to the FS, but are not forbidden 
away from the FS. It is therefore likely that the GNS, 
implied by the global constraint (Eq. (1)) will coincide 
with the FS. 

Generalizations of the above argument to two and 
higher dimensions is straightforward provided the pro- 
jective constraint of no double occupancy in real space 
leads to a finite phase shift for forward scattering on the 
FS and the consequent failure of Fermi liquid theory. In 
particular, the 2D case can also be understood in the 
spirit of Anderson's tomographic Luttinger liquid picture 
H where we have a collection of ID chains in k-space. 

By an exact diagonalisation of finite t— J and Hubbard 
clusters, wc find the gap function Ak by diagonalizing the 
two particle reduced density matrix 

A kk ' = (& k ,_k&k',-k') (4) 

which has the eigenfunction decomposition 



A kk ' = X>«A a (k)A* (k'), (5) 

where A Q and A Q (k) are the a-th eigenvalue and eigen- 
function of the N x N matrix A^w- Here N is the 
number of points in the Brillouin zone and the index 
a orders A Q 's as Ai > A2 > A3 > .... Superconduct- 
ing ODLRO is signalled || by a macroscopic separation 
of the largest eigenvalue Ai from the next one A2, i. e. 
Ai — A2 ~ N. The required gap function Ak is the eigen- 
vector Ai(k) corresponding to the largest eigenvalue Ai. 
For example, for the standard BCS ground state, Ai = N, 
Ai(k) = Mk^k and X a = 0, for a = 2, 3 . . . N. 

We have evaluated A Q and A Q (k) using exact diago- 
nalization of 8-site Hubbard and 16-site t — J chains and 
\/8 x V8 Hubbard and 4 x 4 t — J planes. We chose a 
physically relevant range of ¥ « 5 — 10 and — 0.1 — 0.5. 
For the 16-site t — J chain, we find that Ak has a node 
at k = k s (Fig. 1), and within finite size limitations 
k s = kp, as seen from Fig. 2. In contrast, in a negative- 
U Hubbard model we find no node, as expected (fig. 1). 
In our opinion this is a strong indication for the GNS to 
coincide with the FS in ID. In Figure 3, we display the 
eigenvalue spectrum for the two-particle density matrix 
in order to display the separation of the largest eigen- 
value. For 2D clusters of size v8 x v° and 4x4 with 2 
holes we find that the eigenfunction corresponding to the 
largest eigenvalue, Ai(k), has s-symmetry with a GNS 
close to the FS, as seen from Fig. 4. For the 4x4 case 
with 2 holes, the second eigenvector A2(k), has d x 2_ y 2 
symmetry (Fig. 4). When we look at the eigenvector 
corresponding to the d x i_ y i symmetry, we also find a 
suppressed pairing just outside the FS. This could also 
imply that in the thermodynamic limit, the local con- 
straint may also develop GNS even in the d-state. Note 
that an earlier work of Riera and Young and a recent work 
of Ohta et ai. JIo[| is not inconsistent with the present nu- 
merical results, and in particular, they were not studying 
the behaviour of GNS, which is our main focus. 

We now turn to some of the relevant experimental work 
in this context. From NQR NMR ||,|l) studies, it is clear 
that ^ follows the Korringer law i = kT for T < 10°K, 
in the superconducting state of La2- x Sr x CuO^ for vari- 
ous dopings, suggesting a finite density of quasiparticles 
at the chemical potential. One of the existing explana- 
tions for this is the presence of strong scatterers (close to 
the unitarity limit), which in a d-wave superconducting 
state will lead to gapless excitations fPHl . The number 
of strong scatterers N s needed to explain the required 
density of states at the fermi level is unphysically high 
(n ~ 0.2!). The ubiquitous linear specific heat in the 
superconducting state also seems to survive even in purer 
samples with substantial slope |l2). For La-iCuO^ neu- 
tron scattering also gives an S(q, u>) which indicates the 
presence of gapless excitations in the superconducting 
state. Once again we need a large density of strong 
impurity scatterers to explain the observed density of 
states. Indications for gapless excitations in YBCO is 
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also present in the infrared conductivity studies [[L3| . One 
distinct possibility, consistent with the proposal that we 
make in this letter, is that these gapless excitations are 
intrinsic and are arising from the regions of FS where the 
superconducting gap vanishes on the FS (but are present 
inside and outside). Since there are two parameters in 
the present scenerio, one, the area over which the gap- 
lessness on the FS is there and the other, the slope of the 
gap across the FS on the GNS, we can make a simple fit 
to the experiment, particularly NMR relaxation giving 
a value of about 30% of gapless region on the FS and a 
value of the slope A2 w (when the gap is expressed 
in terms of fermi energy) . As doping increases from zero, 
the FS shrinks dragging part of the GNS with it. (Re- 
call that the GNS coincides with the FS in 2D in our 
numerical result and mean field theory at half- filling) . 
What causes the non-overlap of the GNS with the FS on 
some parts in the BZ? One possible explanation is the 
enhancement of the interlayer pair tunneling matrix ele- 
ment in those directions in k-space. For example, 

as emphasized by Chakravarty et al. Jl5| , tj -.^ is largest 
in the (0, tt) and (tt, 0) directions which enhances pairing 
in those regions of the FS by keeping the GNS away. Our 
suggestion (Fig. 5) incorporates this idea as well as the 
eight nodal points of the recent ARPES data j[6| . What 
we suggest is that the eight points are likely to be nodal 
lines which could be resolved with better experimental 
resolution. (Fedro and Koelling, Chen and Tremblay 
have suggested the possibility of several point nodes (in 
2D) in an extended s-wave state). More importantly, it 
is important to track the GNS whose existence is de- 
manded by a t — J type of modelling. (Finer structure 
in the GNS-FS overlap like possible lumps in the (tt, tt) 
direction could arise due to residual interactions and this 
is incorporated in the figure). 

Figure 6 is an attempt to incorporate d-symmetry in 
our scheme, although at present neither our numerical 
results nor any experiment indicate such a form of the 
gap function. Thermodynamic measurements as well as 
TjT- behavior can hardly distinguish figures 5 and 6. How- 
ever, tunneling measurements can do so and that brings 
us to the issue of the symmetry of the gap. As has been 
pointed out by Stamp and collaborators Q the gap- 
lessness on the FS makes many of the properties simi- 
lar to a gapless situation like a d-wave superconductor. 
Since several tunneling measurements seem to indicate a 
d x 2_ y 2 symmetry, we would like to stress that the gap- 
lessness that we have discussed till now can occur in sys- 
tems with either s or d symmetry. However additional 
point-like nodes (in 2D) as demanded by the symmetry 
will occur in a d-wave situation. 

In conclusion, we have provided some qualitative argu- 
ments with numerical and experimental support to sug- 
gest a partial overlap of the GNS with the FS. It should 
be pointed out that our numerical support for both ID 
and 2D is rather good (though better for ID) and it is 
important to study the shape of the GNS with higher 



experimental accuracy. 
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Figure Captions 

1. Gap function Ai (k) for a 16-site t — J chain at half 
filling (curve A, J = 0.24) and 2 holes (curves B, 
J = 0.08 - 0.32). k F = § , ^ for and 2 holes 
respectively. For comparison, we have also plotted 
the gap function for an 8-site Hubbard chain with 
U=-5 (curve C, 2 holes, k F = \ ). 

2. k s vs. kp for 2, 4, 6 and 8 holes in a 8-site Hubbard 
chain with U = 10, where k s is the node of At- The 
straight line is the k s — kp line as a "guide to the 
eye". 

3. Eigenvalue spectrum of the density matrix for the 
16-site t—J chain. (1) 2 holes, J = 0.24; (2) triplet 
anti-parallel eigenspectrum for 2 holes, J = 0.24. 

4. The topmost two gap functions Ai(k) and A 2 (k) 
for a 4 x 4 plane with two holes, J=0.24. The num- 
bers above and below a square are the coefficients 
of the gap function respectively for the top-most 
(odd-paired) and the next (d-wave) state. 

5. A suggested form of the GNS with s-symmetry. 
The solid line is the FS and the dotted line is the 
GNS. The + and - symbols indicate relative signs 
of the gap function across the GNS. 

6. A possible form of the GNS incorporating d-wave 
symmetry. Conventions are as in Figure 5. 



4 



Figure 1 




Figure 2 
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Figure 3 
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Figure 4 
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Figure 5 
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Figure 6 




